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We examine the entanglement properties of the Yang-Mills theory by calculating a entanglement 
entropy with a ~ 2 using a SU(3) quenched lattice gauge simulation both in the confinement 
and the deconfinement phases. In the confinement phase, the derivative of the a entropy with 
respect to the size / of the subregion, whose entanglement properties are interested in, scales as 
1//^, and a clear discontinuity cannot be found within our statistical errors. The a entropy in 
the deconfinement phase saturates at large /. The saturation value is comparable with the thermal 
entropy of the pure Yang-Mills theory, indicating that the a entropy obeys the volume law at large 
I in the deconfinement phase. 
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1. Introduction 

Various quantum systems show entanglement properties and they receive much attention in 
quantum information theory and condensed matter physics. Entanglement entropy is one of quan- 
tities measuring quantum entanglement. A typical example of the entangled state in quantum me- 
chanical systems is a two spin- 1/2 system in spin singlet state, which is widely used to discuss the 
EPR paradox, one of major topics in quantum physics. Entanglement entropy can be defined not 
only in quantum mechanical systems but also in quantum field theories. 

In quantum mechanical systems, fundamental degrees of freedom are particles and quantum 
entanglement measures how much two or more particles are quantum mechanically correlated with 
each others. In quantum field theories, we focus on quantum entanglement of two or more sibre- 
gions. The entanglement entropy between two subregions, a subregion A of size / and its comple- 
ment B, measures how the spatial subregion in a total system is entangled quantum mechanically 
with its complement. 

Quantum entanglement of ground states has been widely discussed in condensed matter physics 
(for a review, see [1]). For example, the entanglement entropy in the Ising chain model shows a 
divergent behavior at the critical point while it saturates in the non-critical regime. It means that 
the entanglement entropy serves as an order parameter of quantum phase transitions. Therefore, 
the entanglement entropy is a useful quantity to investigate phase structures of quantum systems. 

The entanglement entropy of the pure Yang-Mills theory is particularly interesting. A schematic 
picture of the pure Yang-Mills system is drawn in Fig. 1. The Yang-Mills theory is an asymptot- 
ically free theory and the high energy phenomena in QCD can well be described by gluon and 
quark degrees of freedom using the perturbation theory. At low energies, on the other hand, color 
degrees of freedom are confined in hadrons and the Yang-Mills system is described by the colorless 
hadrons. This may remind us of the deconfinement phase transition; the color degrees of freedom 
are released above the critical temperature, and gluons (and quai^ks in QCD) play a major role as 
effective degrees of freedom while those in the confinement phase are glueballs (or hadrons). Thus, 
one might ask if there is a critical distance scale at which the effective degrees of freedom change 
from colorful objects to colorless objects as the critical temperature of the deconfinement phase 
transition. 

Recently, gauge/gravity correspondence has been extensively studied and it provides a power- 
ful tool to study non-perturbative infrared dynamics of confining gauge theories. Beginning with 
the pioneering work by Ryu and Takayanagi [2], the holographic approach is applied to the calcu- 
lation of the entanglement entropy (for a review on the holographic calculation, see [3]). In this 
approach, the entanglement entropy of gauge theories is obtained by calculating geodesies (mini- 
mal surface bending down to the bulk space) in the gravity side, similar to the calculation of the 
Wilson loop in the holographic approach. The boundary of geodesies coincides the boundary of 
partitioned subsystems in gauge theory side. The entanglement entropy has been studied for vari- 
ous confining backgrounds [4, 5]. It has been argued that the entanglement entropy could exhibit a 
non-analytic behavior with respect to the size / of the subregion; an 0{N^) solution dominates at 
small /, and a /-independent 0(1) solution dominates above some critical length /* (see Fig. 2). This 
indicates that the effective degrees of freedom change from colorful objects to colorless objects, 
and the critical length /* plays a role of the inverse of the critical temperature of the deconfinement 
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Figure 1: Schematic picture of the pure Yang- 
Mill system. 



Figure 2: A typical example of the holograhic 
prediction showing the discontinuity of the en- 
tanglement entropy Sa- 



phase transition. 

The entanglement entropy in SU(2) lattice gauge theory has been studied by Velytsky [6] and 
Bividovich and Polikarpov [7]. In Ref. [6], SU(N) lattice gauge theories are studied in Migdal- 
Kadanoff approximation, and in Ref. [7], SU(2) lattice gauge theory is numerically investigated, 
and there is an indication that the derivative of the entanglement entropy shows a discontinuous 
change at some critical length scale /* and it vanishes at large /. 

In this study, we investigate a entanglement entropy in SU(3) pure Yang-Mills theory using 
lattice Monte Carlo simulations. Instead of directly calculating the entropy, we adopt numerical 
technique to evaluate the entanglement entropy, which has also been used in [7] (originally pro- 
posed in [8, 9] in order to calculate the pressure in the deconfined phase). 



Entanglement entropy of a pure state is defined 
as follows. We divide the total system into subregion A 
and its complement B (see Fig. 3). Let / be the size of 
the system A in the x direction. The density matrix of the 
system is p = At zero temperature, the ground 

state is a pure state and the von Neumann entropy of the 
system is zero. The reduced density matrix obtained by 
tracing out the degrees of freedom in the region B, 

Figure 3: The complementary regions = Trg p = Trg , (2. 1) 

A and B separated by an imaginary 

boundary at x = /. describes the density matrix for an observer who can only 

access to the subregion A. Although we start off with a 
pure state with vanishing von Neumann entropy, the state corresponding to the reduced density 
matrix is generally a mixed state, p^ contains the information on the quantum degrees of freedom 
traced out. The entanglement entropy is defined as the von Neumann entropy of the reduced density 




3 



Entanglement entropy 



Y. Nakagawa 



matrix, 

Sa = -TrpAlnpA. 
Some properties of the entanglement entropy can be found in [10]. 



(2.2) 



3. Replica method 



Z(/,a) = 




Figure 4: Schematic picture 
for the system with a cuts 
in X — f plane. In the region 
A (B), the periodic boundary 
condition is imposed with 
the period a/3 (j3). 



In order to evaluate the entanglement entropy, we apply the 
replica trick. The detail of the derivation is given in [11]. The point 
is that the entanglement entropy defined in Eq. (2.2) can be repre- 
sented in the form, Sa = — limce^.i d/dalnTvAP^. The trace of the 
a-th power of the reduced density matrix Pa is given by the ratio of 
the partition functions. 



^ „ Z(l,a) 



(3.1) 



Here Z(/,a) is the partition function of the system having special 
topology, the a-sheeted Riemann surface, and Z = Z{a = 1). The 
field variables in the region A is periodically identified with the inter- 
val aj3 (j8 = 1/r is the lattice size in the temporal direction) while 
in the region B the periodic boundary condition is imposed with the 
period jS (see Fig. 4). 

The entanglement entropy is then given by 



SAil) 



d 

lim ^- — In 



Z{l,a) 



(3.2) 



The derivative of Sa{1) with respect to /, which is free of the ultraviolet divergence, can be expressed 
as follows; 



dl 



d 
dl 



lim In 

a^l da 



Z{l,a] 



lim ^^F[l,a]. 
«->! ol da 



(3.3) 



That is, in order to calculate dSA/dl, we evaluate the free energy of the system having a cuts as 
is depicted in Fig. 4, take the derivative with respect to a and /, and then take the limit a — 1. 
Thus, the evaluation of the entanglement entropy is reduced to the calculation the free energy of 
the system with a cuts. 



4. Lattice setup and observables 

In the lattice simulations, the derivative in Eq. (3.3) is replaced by the finite difference, and we 
estimate the derivative by 

hm — —F A, a ^ ^ lim f /,« + 1 -f /,« ^ ^ ^- ^. (4.1) 

«->! dl da dl a 

We note that dF[l,a = \\/dl drops out since F[l,a = 1] does not depend on /. The difference 
of the free energies can be evaluated numerically by introducing an 'interpolating action' which 
interpolates two actions corresponding to two free energies [8, 9], Sint = (1 — Y)Si[U] + ySi+aiU]. 
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Si and Si+a represent the actions corresponding to F[l, a = 2] and F[l + a, a = 2] in Eq. (4. 1). It is 
easy to show that 

F[l + a,a = 2]-F[l,a = 2] = - j\y^\nZ{l,y) = j\y{Si+a[U]-S,[U])^. (4.2) 

Here (•)^ refers to the Monte Carlo average with the interpolating action Smi- Therefore, the a = 2 
entanglement entropy can be evaluated numerically by updating gauge configurations with the 
interpolating action on the lattice with a = 2 cuts, calculating the action differences for various y, 
and performing a numerical integration over y. In order to evaluate the integral in Eq. (4.2), we 
calculated the action differences from 7 = to 1 by the step 0. 1, and employed the Simpson's rule 
to evaluate the integration numerically, which interpolates neighboring points by a quadratic curve. 

We adopt the heat-bath Monte Carlo technique with the standard plaquette action to generate 
lattice configurations. First 5000 sweeps are discarded for thermalization, and the measurement 
has been done every 100 sweeps. The number of configurations for each jS and lattice size is about 
3000 to 8000. The statistical errors are estimated by the jackknife method. 



5. Simulation results 



5.1 a entanglement entropy at zero temperature 

The derivative of 5a (/) with respect to / in the confinement phase is plotted in Fig. 5. dSA{l)/dl 
is normalized by the area of the common boundary of the two subregions, \dA\. We observe that 
data on 12^ and 16^ agree within statistical errors. This implies that the derivative of the a = 2 
entanglement entropy is proportional to the area of the boundary, namely, entanglement entropy 
obeys an area law at zero temperature. 

In the small / region, the a = 2 entanglement entropy is expected to scale as 1 //^ from the 
dimensional analysis. That is, dSA/dl behaves as \/P at small /. This behavior is exactly what the 
entanglement entropy in conformal field theory in (3+l)-dimensional spacetime shows. In order to 
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Figure 5: (1 /\dA\)dSA/dl in the confinement phase. The dashed curve is the fit of the data by the function 
c//'' with the fitted values c = 0. 149(48), c/ = 3.06(20). The right panel shows the zoom up of the left panel 
to make near-zero region more visible. 
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confirm this, we fitted data with the function BSa/SI = c{\/lY , and we obtain c = 0. 149(48), (i = 
3.06(20), /ndf = 0.192. The fitted function is plotted in Fig. 5 by the dashed curve. Our result 
does not show a clear sign of the discontinuity in dSA/dl. Since the derivative of Sa rapidly 
decreases, the signal-to-noise ratio becomes quite small at large / and it is very difficult to locate 
the critical length of the entanglement entropy numerically, even if it exists. It can be safely stated 
that our results exclude the possibility of the existence of the critical length at below 0.4 [fm]. 



5.2 a entanglement entropy below and above the critical temperature 

The left panel of Fig. 6 shows the derivative of the a entanglement entropy below the critical 
temperature. The fitted function at zero temperature, {\/\dA\)dSA/dl = 0.149//^ *'^, is drawn by 
the dashed curve. We observe that the data agree with the fitted function of the zero temperature re- 
sult. This indicates that the a entanglement entropy does not show a clear temperature dependence 
below the critical temperature. 

The numerical result above the critical temperature is given in the right panel of Fig. 6. We see 
that dSA/dl does not approach zero but saturates at large /. At zero temperature, the ground state is 
a pure state and the von Neumann entropy is zero. By contrast, the ground state at finite temperature 
is a mixed state and the von Neumann entropy (thermal entropy) takes a finite value. This means 
that at finite temperature, the entanglement entropy measures not only the quantum mechanical 
con^elation between the two spatial subregions but also the thermal entropy of the subregion. Since 
the thermal entropy of the SU(3) Yang-Mills theory rapidly increases in the vicinity of the critical 
temperature, the saturation value of the entanglement entropy may be considered as the thermal 
entropy of the subregion A. We note that the asymptotic behavior of {\/\dA\)dSA/dl implies that 
the entanglement entropy obeys the volume law at lai^ge / above the critical temperature. 

We fitted the data with the function a/P + b, and we obtain 

= 0.180(16), Z7= 14.1(8) (r/r,~1.44) 
= 0.187(14), Z7 = 61.8(15) (r/r,~2.02). 
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6: Left panel: (1 /\dA\)dSA/ dl below the critical temperature. The dashed curve is the fitted function 
temperature. Right panel: {\ /\dA\)dSA/ dl above the critical temperature. The dotted curves show 
of the data by the function c/fi +d. 
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We note that the coefficient of the l/P term agrees with each other. In order to compare the 
asymptotic value of dSA/dl to the thermal entropy, we estimated the thermal entropy SLtT /Tc ^ 
1.44 and 2.02 by reading the values off the figure in Ref. [12]. A rough estimate gives 

.= 17 (r/r,~i.44) 

s = 56 {T/Tc ^2.02). 
These are comparable with the asymptotic values of dSA/dl. 

6. Summary and conclusion 

We studied the a entanglement entropy of the Yang-Mills vacuum with a = 2 using lattice 
Monte Carlo simulations. The entanglement entropy measures the quantum correlation between 
spatial subregions. We find that at zero temperature the derivative of the a = 2 entropy with re- 
spect to / is well fitted by the function c/l'' with d=3.06(20). The exponent d is consistent with 
that in the conformal field theory. A clear discontinuity in dSA/dl was not observed within the sta- 
tistical errors, which is arguded in the models of the gauge/gravity correspondence. Furthermore, 
we observe that the a entropy is almost temperature independent below the critical temperature. 
Above the critical temperature, a = 2 entropy does not approach zero but saturates at large I. Since 
the ground state of the finite temperature system is a mixed state, this implies that the entangle- 
ment entropy measures not only the coiTclation between the spatial subregions but also the thermal 
entropy of the subregion, which dominates at large /. Indeed, our fitted result of the asymptotic 
values is comparable with the thermal entropy of the pure SU(3) Yang-Mills theory. 
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